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Abstract 

We compute the S'-matrix of the Tricritical Ising Model perturbed by the sub- 
leading magnetic operator using Smirnov's RSOS reduction of the Izergin-Korepin 
model. The massive model contains kink excitations which interpolate between 
two degenerate asymmetric vacua. As a consequence of the different structure of 
the two vacua, the crossing symmetry is implemented in a non-trivial way. We use 
finite-size techniques to compare our results with the numerical data obtained by 
the Truncated Conformal Space Approach and find good agreement. 
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1 Introduction 



It has been pointed out by Zamolodchikov that certain deformations of minimal models 
of Conformal Field Theories (CFT) produce integrable massive field theories, which are 
characterized on mass-shell by a factorizable S'-matrix [|T|. A particularly interesting 
situation occurs for the "Tricritical Ising Model" (TIM) at the fixed point perturbed 
by the operator $2,1. This field has anomalous dimensions (A, A) = (j^, ^) and is 
identified with the sub-leading magnetic operator[|, odd with respect to the Z2 spin- 
reversal transformation. Hence, this deformation explicitly breaks the Z2 symmetry of 
the tricritical point and the corresponding massive theory can exhibit "^■^-property" , 
i.e. the absence of a conserved current of spin 3 and the possibility to form a bound 
state through the process AA A AA. The counting argument supports this 
picture, giving for the spin of the conserved currents the values s=(l,5,7,ll,13) 0]. The 
interesting features of this massive field theory have been first outlined in ref. [Q] where 
the model was studied by the "Truncated Conformal Space Approach" (TCSA), proposed 
in ref. [0]. This approach consists in the diagonalization of the perturbed Hamiltonian 

E = EcFT^\ J <^2,i{x)dx (1.1) 

on a stripQ of width R, truncated at a certain level of the Hilbert space defined by the 
conformal field theory at the fixed point. The lowest energy levels are given in fig. 1. The 
theory presents two degenerate ground states (which correspond to the minima of the 
asymmetric double-well Landau-Ginzburg potential in fig. 2) and a single excitation B of 
mass m below the threshold at 2m. The double degeneracy of the vacuum permits two 
fundamental kink configurations | K^) and | K^) and, possibly, bound states thereof. If 
the two vacua were related by a symmetry transformation, one would expect a double 
degeneracy of the breather-like bound state | B) in the infrared regime R 00. However, 
the absence of a Z2 symmetry makes it possible that in this case only one of the two 

asymptotic states | K^KJ) or | K^K^) is coupled to the bound state | B). 

^In the Landau-Ginzburg theory, the TIM represents the universahty class of the </j^ model. The 
operator $2,1 is represented by the third-order field (^^. 

^We consider in the following only the case of periodic boundary condition. 



A scattering theory for such model with asymptotic states | | KJ) and | B) 

was first conjectured in 0]. On the other hand, a general framework for the $1.2 and 
$2,1 deformations of CFT has been recently proposed in p. It is based on the RSOS 
reduction of the Izergin-Korepin model 0. In section 3 we explicitly work out the RSOS- 
like S-matrix in the case of the TIM perturbed by the subleading magnetic operator. The 
scattering theory we obtain is different from the one of ref. although both of them 
have the following features: 

1. the existence of two fundamental kink configurations with mass m; 

2. the appearance of only one bound state of the above, with the same mass m. 

Hence, both scattering theories give rise to a picture that qualitatively agrees with the 
results of ref. 0] . Therefore a more detailed analysis is required in order to decide which 
of the two theories is appropriate for the description of the scaling region of the TIM 
in presence of a third-order magnetic perturbation. As suggested in ref. 0|, possible 
insight comes from the study of finite-size effects for the theory defined on the cylinder 
1^. These finite-size effects, which can be directly related to the scattering data, control 
the exponential decay of the one-particle energy level to its asymptotic value. 

The paper is organized as follows: in sect. 2 we outline Smirnov's RSOS reduction of 
the Izergin-Korepin model. In sect. 3 we present our results on the S'-matrix of the $2,1 
perturbation of TIM. Sect. 4 contains the analysis of the phase shifts and the asymptotic 
behaviour of the RSOS S'-matrix. Sect. 5 consists of a short summary of Zamolodchikov's 
proposal for the S- matrix |^. This proposal, together with ours, will be checked against 
numerical data in sect. 6, where we perform an accurate analysis of the energy levels 
obtained from the TCSA using the finite-size theory. Our conclusions are in sect. 7. 

2 Smirnov's RSOS reduction of the Izergin-Korepin 
model 

In a recent paper, Smirnov has related the $1^2 and $2,1 deformations of a CFT to a RSOS 
projection of the i?-matrix of the Izergin-Korepin (IK) model |^. His results can be 



summarized as follows. The i?-matrix of the IK model cannot be directly interpreted 
as the S'-matrix of the massive excitations of the perturbed CFT because, introducing a 
Hilbert space structure, it does not satisfy the unitarity requirement. But it is possible 
to use the quantum SL{2)q symmetry of the IK model to study the RSOS restriction 
of the Hilbert space. This happens when g is a pth root of unity. The RSOS reduction 
preserves the locality of an invariant set of operators and yields S'-matrices which have a 
sensible physical interpretation. The RSOS states appearing in the reduced model 

I ji, fci, I ai I /32, j2, A;2, . . . I a^-i I Pn,jn,kn) (2.1) 

are characterized by their rapidity (3i, by their type k (which distinguish the kinks from 
the breathers), by their SL{2)q spin j and by a string of integer numbers aj, satisfying 

p — 2 

Qi < — ^ , I ctfc - 1 |< Ofc+i < min(afc + - 3 - a^) . (2.2) 

In the case of $2,1 deformation, using the following parameterization of the central charge 
of the original CFT 

c=l-&[^ + ^-2j , (2.3) 

the g-parameter is given by 

q = exp(2i7) . (2.4) 

For the unitary minimal CFT, 7 = vr ^qj^ (r = 3,4, . . .). It is also convenient to define 
the quantity 

The S-matrix of the RSOS states is given by p[ 



SoiPk — Pk+l) 




(2.6) 



- fexp l-^iPk+i - Pk] 



A ^-(Cafe + l+Cafc_l-Cafc-C„J^+3)/2 



Herein, Ca are the Casimir of the representation a, Ca = a{a + 1) and the expression of 
the 6j-symbols is that in ref. [0. So{l3) has the following integral representation 



sinh-(/?-m)sinh-(^/3- — 



X exp —2i 



dx sin Px sinh ^ cosh (| - |) 



X 



X 



cosh ^ sinh ^ 



(2.7) 



3 RSOS S'-matrix of <|)2,i perturbed TIM 

In the TIM perturbed by the sub leading magnetization operator, r = 4 and ^ = 



From eq. ( |2.2| ), the only possible values of are and 1 and the one-particle states are 
the vectors: | Kqi), \ Kiq) and | Ku). All of them have the same mass m. Notice that 
the state | Kqq) is not allowed. A basis for the two-particle asymptotic states is 



liToi^io), \KnKn), | i^iii^io), | i^io^^oi) • 

The scattering processes are 



(3.1) 
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(3.2) 




= sllip^ 


-P2) 


|i^ll(/52)i^ll(/9l))- 


h sllip. 


-P2) 


1 K,,{h 


)i^oi(/3i)) 




= s^APi 


-P2) 


\K,,{(32)Koi{l3i))- 




-P2) 


1 ^11 (/?2, 


)^ii(/5i)) 



Explicitly, the above amplitudes are given by 



9 ^ .TT^ 



(3.3.a) 



9 . .TT^ 



oN<Ql = S'M = -^5o(/?) sinh(-/3 + z^) 



(3.3.b) 



\-^/ ,1 i sin(f) /g 



ill 



^9^ .27r 



(3.3.C) 



iXi 






s'M = - 



sm I ^ 



sinh Q/?) (3.3.d) 



/ \ 



?■ sin (I) /Q 27r 



(3.3.e) 



The function 5'o(/9) can actually be computed. It is given by 



So(P) = -(sinh^(;3-ix)sinh^(/J-^) 

x»(A-i)™(A+i)»(Al) 
x<(ft^)t(ft-^)*(Al)*(ft-i), 



where 



w(/5,a;) 



mx 



sinh + 

sinh — iTTx) 

sinh + i'Kx) 

sinh ^(/S — iTTx) 



(3.4) 



It is easy to check the unitarity equations: 



S',li-P) + SlliP) Slli-P) = 1 ; 
Sl'A^) + Sim 5n°(-/5) = ; (3.5) 

Sim simp) = i- 

^oo(/5)^ooM) = l- 

An interesting property of this S matrices is that the crossing symmetry occurs in a 
non-trivial way, i.e. 

Sllitn-P) = SlliP); 

5°°(z7r-/?) = a' Sim-, (3.6) 
S°|(m-/?) = aS'M; 



where 



(3.7) 



Af). 

and s{x) = sin(7rx). 

The above crossing-symmetry relations may be seen as due to a non-trivial charge 
conjugation operator (see also 0). In most cases, the charge conjugation is implemented 
trivially, i.e. with a = ±1 in eq. ( [3.61 ). Here, the asymmetric Landau-Ginzburg potential 
distinguishes between the two vacua and gives rise to the value 

The amplitudes (3.3) are periodic along the imaginary axis of f3 with period 10 iri. The 
whole structure of poles and zeros is quite rich. On the physical sheet, < Im (3 < in, 
the poles of the ^-matrix are located at /? = ^ and /5 = y (fig. 3). The first pole 
corresponds to a bound state in the direct channel while the second one is the singularity 
due to the particle exchanged in the crossed process. The residues at /5 = ^ are given 
by 

11/ 



2 



ri = Res^^2|i Sqq{(3) = ; 

r2 = Res^^2p S^l{(3) = i (^j|fjj ^5 
ra = Res^^22i Sll{P)=iuj; (3.8) 



where 




Their numerical values are collected in Table 1. 

In the amplitude Sqq there is no bound state in the direct channel but only the 
singularity coming from to the state | Ku) exchanged in the t-channel. This is easily 
seen from Fig. 4 where we stretch the original amplitudes along the vertical direction (s- 
channel) and along the horizontal one (t-channel). Since the state | Kqq) is not physical, 
the residue in the direct channel is zero. In the amplitude SqI we have the bound state 
I Kqi) in the direct channel and the singularity due to | Ku) in the crossed channel. In 
5*^^, the state | Ku) appears as a bound state in both channels. In Sil the situation is 
reversed with respect to that of S^l, as it should be from the crossing symmetry property 
( ^.6| ): the state | Ku) appears in the t-channel and | Kqi) in the direct channel. Finally, 
in there is the bound state | Ku) in the direct channel but the residue on the t- 
channel pole is zero, again because | Kqq) is unphysical. This situation is, of course, that 
obtained by applying crossing to Sq^. 



4 Energy levels, phase shifts and generalized statis- 
tics 

The one-particle line a of fig. (l.a) corresponds to the state | Ku)- This energy level 
is not doubly degenerate because the state | Kqq) is forbidden by the RSOS selection 
rules, eq. (|2.2|). With periodic boundary conditions, the kink states | Kqi) and | i^io) are 
projected out and | Ku) is the only one-particle state that can appear in the spectrum. 

Consider the threshold line. On a strip with periodic boundary conditions we ex- 
pect this energy level to be doubly degenerate. This because the states | KqiKiq) and 



I KiqKqi) are identified under such boundary conditions whereas the state | KuKu) re- 
mains distinct. The conformal operators creating these states in the u.v. hmit are $x x 

16 '16 

and $A 6_, respectively (see hues marked b and c in fig. l.a). We have checked that for 
large values of R these two lines approach each other faster than 1/r (as would be the 
case if the lowest line b was the only threshold line and c a line of momentum). However 
the truncation effects already present in this region prevents us from showing that they 
really approach each other exponentially^. In the following we assume that ideally these 
two levels are exponentially degenerate in the infrared limit. 

The identification of the threshold lines described above holds only for the static 
configuration of two kinks with zero relative momentum. The situation is indeed different 
for the lines of momentum which approach the threshold. In fact, the ^-matrix acting 
on I KoiKio) can only produce | KqiKiq) as final state whereas | KiqKqi) and | KnKn) 
can mix through the processes of eq. (3]^). In order to determine the pattern of the 
energy levels obtained from TCSA and to relate the scattering processes to the data of 
the original unperturbed CFT (along the line suggested in |jl2[); we would need a higher- 
level Bethe ansatz technique. This is because our actual situation deals with kink-like 
excitations in contrast to that of ref. [jl2| which considers only diagonal, breather-like S 
matrices. The Bethe-ansatz technique gets quite complicated in the case of a S'-matrix 



with kink excitations. For the moment, it has been applied in few cases |]T3|, In 
the light of these difficulties, we prefer here not to pursue such a program and instead 
concentrate on some properties of the phase shifts and generalized statistics which arises 
for kink excitations. 

For real values of P, the amplitudes Sqq{P) and SqKP) are numbers of modulus 1. It 
^Consider, for instance, figure 10 of ref. |^ (the case of low-temperature phase of thermal perturbation 
of the TIM) . One observes that the onset of the exponential approach of degenerate excited levels usually 
occurs quite far from the value of R at which the two ground state energies coincide. As we discuss in 
sect 6, the situation of the subleading magnetic deformation of the TIM is even worse from a numerical 
point of view, because the anomalous dimension of the subleading magnetic operator is approximately 
1/2. In this case, it is therefore possible that the onset of the exponential approach of the higher levels 
takes place in a region of R strongly dominated by truncation effects. 



is therefore convenient to define the following phase shifts 



S'M ^ e2^^«(^); (4.1) 



The non-diagonal sector of the scattering processes is characterized by the 2x2 symmetric 
^-matrix 

' Sim Sim 
Sim Sim 



(4.2) 

We can define the corresponding phase shifts by diagonalizing the matrix ([4.21). The 



eigenvalues turn out to be the same functions in ([4.1|) , 

/ 



g2j5o(/3) 






(4.3) 



The phase shifts, for positive values of /3, are shown in fig. 5. Asymptotically, they have 
the following limits 

hm e^^'^^^^ = e±¥; (4.4) 
hm e2^^^(^) = e±¥. 

There is a striking difference between the two phase shifts: while Si{(3) is a monotonic 
decreasing function, starting from its value at zero energy 5i(0) = |, 5o(/5) shows a 
maximum for ~ f and then decreases to its asymptotic value Its values are always 
larger that 5o(0) = f . Such different behaviour of the phase shifts is related to the 
presence of a zero very close to the real axis in the amplitude e^*'^"''^\ i.e. at /3 = i^. 
This zero competes with the pole at /3 = i| in creating a maximum in the phase shift. 
Similar behaviour also occurs in non-relativistic cases [rU| and in the case of breather-like 



5" matrices which contains zeros [|rT|]. The presence of such a zero is deeply related to the 



absence of the pole in the s-channel of the amplitude e^^^°^^\ For the amplitude e'^^^'^^l^\ 
the zero is located at (3 = ^ (between the two poles) and therefore its contribution to 
the phase shift is damped with respect to that one coming from the poles. The net result 
is a monotonic decreasing phase shift. 



Coming back to the 2x2 S- matrix of eq. ( |4.2| ), a basis of eigenvectors is given by 
I 0i(/5i)0i(/?2)) = A(/3i2) (I KumK,,{P2)) + Xi{Pi2) I K.omKoim)) ; (4.5) 

where A{i3i2) is a normahzation factor. In the asymptotic regime f3 —>■ 00 

27vi / o Gtti \ 

e 5 a + e 5 

XI = ^ (4.6) 

a 

X2 = ^ , 

a 

and the probabihty Piooi to find a state | KiqKqi) in the vector | 0202) w.r.t. the 
probabihty Puu to find a state | KuKu) is given by the golden ratio 

^1001 1 /7r\ 

-7 = 2 cos - . (4.7) 



Pun a2 V5 

For the state | we have 



a'= . (4. 



nil 



2cos(f) 



The "kinks" 0i and 02 have the generahzed bihnear commutation relation [|T5| 0, [17 



0,(t, x)0,(t, y) = 0,(t, 2/)0i(t, x) e2-^-^(^-^) . (4.9) 

The generalized "spin" is a parameter related to the asymptotic behaviour of the 
S-matrix. A consistent assignment is given by 

_ 3 _ 5o(oo) 
•Sii — - — ; 

TT 

S12 = ; (4.10) 



S22 = T7T 



3 5i(oo) 



10 7r 

The implications of these generalized statistics will be discussed elsewhere. Here we only 
notice that, interesting enough, the previous monodromy properties are those of the chiral 
field \I' = $iL of the original CFT of the TIM. The operator product expansion of \l/ 
with itself reads 

^iz)^iO) = 4 1 + ^(0) + . . . (4.11) 

Z5 Z5 



where Cq,^^^^ is the structure constant of the OPE algebra. Moving z around the origin, 
z — >■ e'^'^'^z, the phase acquired from the first term on the right hand side of ([4.11|) 



comes 



from the conformal dimension of the operator \1/ itself. In contrast, the phase obtained 
from the second term is due to the insertion of an additional operator \E'. A similar 
structure appears in the scattering processes of the "kinks" (pf. in the amplitude of the 
kink 01 there is no bound state in the s-channel (corresponding to the "identity term" in 
( [4.11| )) whereas in the amplitude of 02 a kink can be created as a bound state for P = ^ 
(corresponding to the "\E' term" in ( [4.11|) ). In the ultraviolet limit, the fields 0i should 
give rise to the operator "^{z), similarly to the case analyzed in The actual proof 
requires the analysis of the form factors and will be given elsewhere. 



5 Zamolodchikov's S'-matrix for ^2,1 perturbed TIM 

The problem of finding a theoretical explanation for the energy levels of the $2,1 perturbed 
TIM was first discussed in a remarkable paper by Zamolodchikov [Q. In his notation, the 
one particle states are given by 

\K+), \K^), \B), (5.1) 

which we can identify with our | Kiq), \ Kqi) and | Ku), respectively. The two-particle 
amplitudes of the scattering processes were defined in to be 



1 B{(3,)B{P,)) 


= a{Pi 


-P2) 


1 B{h)B{Pi)) + 6(A - /32) 1 K^{P2)K.{Pi)) 




= c(/?i 


-^2) 


1 i^_(/52)5(/?l)) 


1 B{(3,)K+{(32)) 


= c(/3i 




1 B{f32)K+{(3,)) 




= d{(3. 


-P2) 


1 ir_(/32)i^+(/5i)) 


K^{(3,)K_{(32)) 


= e(A 


-P2) 


1 K+(/32)i^-(/3i)) + KPi - P2) 1 B{P2)B{f3i)) 



They are in correspondence with those of eq. (|3.2| ) if we make the following assign- 
ments 



m 

c{f3) 
eif3) 



(5.3) 



In order to solve the Yang-Baxter equations which ensure the factorization of the scat- 
tering processes, Zamolodchikov noticed that the above amphtudes coincide with the 
definitions of the Boltzmann weights of the "Hard Square Lattice Gas". Therefore, he 
borrowed Baxter's solution |Jl8|] in the case where it reduces to trigonometric form 

(f + A/?) 



a(/5) 
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c{(3) 
e{(3) 



sm 



R{/3) 



-513 



(f ) 

sin (A/?) 

^sin (f ) (f ) 
in (f - A/3 



-R{(3) 



-25/3 



sm 

sm 
sin 



-(f) 
sin (f ) 



R{P) ■ 
■R{(3) ■ 



(5.4) 



sin - A/3) 



sm 



Here 5 and A are arbitrary parameters and R{I3) is an arbitrary function. In order to fix 
completely the amplitudes, Zamolodchikov imposed the following requirements: 

1. the unitarity conditions, eqs. (|3.5| ); 

2. the absence of a pole in the direct channel of the amplitude d{[3); 

3. crossing symmetry, implemented in the following form 

a(/5) = a{ir[ — (3) ; 

6(/?) = c(m-/5); (5.5) 
d{p) = e{m - 13) . 



The final form of tlie S matrices is given by 

,.,,sin(2^fS)si„(-±«) 



a{/3) = e 



sin (^) sin (2^±6i^) ' 
sin sin (2^) 



sm I I sm 



sin (^) sin (2zi±6l^) 
"^^^ ^ " sin (^) sin ' ^^-^^ 



e(/3) = e- 



-2<5/3 



Sin (3^±6i^) 



sm 
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-25(i7r-/3) 



sm 



Herein 5 is an imaginary number satisfying 



(!) 

All amplitudes but have a simple pole at /3 = Their residues are given by 



5, 



' (D) 



5(^(1 



13 

5y 



' (I)) 



T3 = Res^^2^ c{(3) = i - ^ ) ^^ ; (5.8) 
T4 = Res^^27!i = ; 

r5 = Res^^22H e(/3) = i | ^ . 

Their numerical values are collected in Table 2. 

In the asymptotic limit /3 — > 00, all amplitudes but a(/3) have an oscillating behaviour 



c(/3) ~ e-'^e-'-^ ■ (5.9) 

Such oscillating behaviour was also found by the same author for a scattering model with 
Z4 symmetry |19|. There, it was suggested that in the ultraviolet limit, the Z4 theory 



has a limit cycle. Whether or not this is also the case for the S-matrix proposed by 
Zamolodchikov for the subleading magnetic perturbation of TIM, it is not clear to us 
that it is possible to match such oscillating behaviour to either a definite CFT in the 
ultraviolet limit, or a choice of generalized statistics for the kinks. 

6 Finite-size effects 

The study of the scaling region around a fixed point is simplified along those directions 
(in the space of coupling constants) which define an integrable massive field theory. In 
such cases, the knowledge of the theory on mass-shell (the S'-matrix) makes it possible 
to characterize completely the dynamics even off-shell. In particular, using the Thermo- 
dynamical Bethe Ansatz method |^| one could compute the ground-state energy Eq{R) 



for the theory on a cylinder of width R. If this computation could be performed for 
the subleading magnetic perturbation of TIM, it would become easy to decide which of 
the two proposed scattering theories is the correct one. Unfortunately, the Bethe-ansatz 
technique has not been extented to the case of an S'-matrix with kink excitations. Only 
few examples have been worked out |TB|, |T^. However, we can get around this difficulty 
using the "Truncated Conformal Space Approach" (TCSA) p|, |, |22| and the predictions 
of finite-size theory |]. 

The TCSA allows us to study the crossover from massless to massive behaviour in 
a theory with the space coordinate compactified on a circle of radius R. The method 
consists in truncating the infinite-dimensional Hilbert space of the CFT up to a level A 
in the Verma modules. Then, the off-critical Hamiltonian 



H{\,R) = Ho{R) + V{R) 



(6.1) 



is numerically diagonalized. Here, Ho{R) is the Hamiltonian of the fixed point ^1 



HoiR) = ^ (Lo + Lo-^) , (6.2) 



i? V ^ " 12 
and V{R) is the interaction term given by the perturbation 



V{R) = A f'^^rsix) dx . (6.3) 
Jo 

The matrix elements of V{R) are computed in terms of the three-point functions of the 
scaling fields of the fixed point. An efficient algorithm has been developed for performing 
such computation [^]. In our case, the truncation is fixed at level 5 in the Verma 
modules. The parameter A in ( |6.3|) is a dimensionful coupling constant, related to the 
mass scale of the perturbed theory 

[A] = m2-2^--= . (6.4) 
In the following we fix the mass scale by A = 1. 



The energy levels Ei[R) of the Hamiltonian ( |6.1|) have the scaling form 



Ei{R) = ^Fi{mR) , (6.5) 



with the asymptotic behaviour 



tQirn^R + Mi , mR > 1 . 

Here, Aj is the anomalous dimension of some scaling field in the ultraviolet regime and 
Mi is the (multi)particle mass term in the infrared limit. However, the above infrared 
asymptotic behaviour holds only in the ideal situation when the truncation parameter 
A goes to infinity. In practice, for finite A, the linear behaviour of eq. (|6.6| ) is realized 
only within a finite region of the R axis. The large R behaviour is dictated by truncation 
effects. In order to find the physical regions, we make use of the following parameter 
(introduced in ref. [0]) 

d log Ei{R) 
d log it 

The parameter pi is between the values pi = —1 (in the ultraviolet region) and pi = 1 (in 
the infrared one). In the limit of large R (the truncation-dominated regime), pi = 1 — 2 A. 



The "window" in R where the hnear infrared behaviour holds depends upon the 
perturbing field and, for the case of operators with anomalous dimension A > i, it can 
be completely shrunk away. This phenomenon is related to the divergences which appear 
in a perturbative expansion of the Hamiltonian ( |6.1| ), which must be renormalized. Under 
these circumstances, it is more convenient to consider the differences of energies, which 
are not renormalized. 

In the case of the subleading magnetic perturbation of TIM, the anomalous dimension 
of $2,1 is A = ^, which is near |. Looking at fig. (l.a), we see that the onset of the 
infrared region of the two lowest levels is around i? ~ 2 and persists only for few units 
in R. In this region one can check that they approach each other exponentially 

Ei-Eq^ e"""^ , (6.8) 

and extract in this way the numerical mass of the kinks 

m = 0.98 ±0.02. (6.9) 

From fig. 6, we see that for the third level, that of one-particle state, the ultraviolet 
behaviour extends till R ~ 0.5. The crossover region is in the interval 0.5 < R < 2. 
Beyond this interval, the infrared regime begins but the "window" of infrared behaviour 
is quite narrow and is in the vicinity of i? ~ 3. In such small region, it becomes hard 
to extract any sensible result. To overcome this difficulty, it is better to consider the 
differences of energies with respect to those of the degenerate ground states. From fig. 
(l.b) one can easily read off the mass-gap and see that is consistent with the value 
extracted from the exponential approach of the two lowest levels. In fig. (l.b), the third 
line defines the threshold, with a mass-gap 2m. 

In the ideal situation A oo, the crossover between the intermediate region [mR ~ 
1) and the infrared one {mR ^ 1) is controlled by off-mass shell effects and has an 
exponential behaviour. The computation of these finite-size corrections has been put 
forward by Liischer [||. In the case of one-particle state, there are two leading off- 
mass-shell contributions coming from the processes shown in fig. 7. The first correction 
involves the on-mass-shell three-particle vertex F, which is extracted from the residue at 



/5 = 2|i of the amplitudes Sll{(3) (in the case of RSOS S-matrix) and a{f3) (in the case 
of Zamolodchikov's S'-matrix). The second correction comes from an integral over the 
momentum of the intermediate virtual particle, interacting via the S'-matrix {Sll{f3) for 
the RSOS S'-matrix and a(/5) for that of Zamolodchikov). Therefore we have 

AE{R) ^ E^iR) - Eo{R) = m + z exp (^-^^^^ (6.10) 



-m 



We have done the following. First we have computed numerically the integral on the 
intermediate particles in both cases of RSOS and Zamolodchikov's S-matrix and we have 
subtracted it from the numerical data obtained from the TCSA. After this subtraction, 
we have made a fit of the data with a function of the form 

G{R) =A + Be-^"^^ + Ce"'"^ . (6.11) 

The first term should correspond to the mass term. The coefficient of second one is the 
quantity we need in order to extract the residue of the S'-matrix at /5 = ^ 




B = iRes^=2^ < . (6.12) 



The third term is a subleading correction related to the asymptotic approach of the lowest 
levels of our TCSA data to the theoretical vacuum energy Eq{R). 

In the case of RSOS S'-matrix, the best fit gives the following values 

A = 0.97 ±0.02; 

B = -0.29 ±0.02; (6.13) 
C = -0.36 ±0.02. 

The corresponding curve is drawn in fig. 8, together with the data obtained from TCSA. 
The mass term agrees with our previous calculation (eq. ( |6.9| )). The second term gives 
for the residue at /3 = ^ the value 0.34 ±0.02. This is consistent with that of the RSOS 
S-matrix. In our fit procedure, the value of the residue we extracted through (|6.12|) is 
stable with respect to small variation of the mass value. Increasing (decreasing) m, B 



increases (decreases) as well, in such a way that the residue takes the same value (into 
the numerical errors). This a pleasant situation because it permits an iterative procedure 
for finding the best fit of the data: one can start with a trial value for m (let's say m = 1) 
and plug it into ( |6.11| ). From the A-term which comes out from the fit, one gets a new 
determination of the mass m that can be again inserted into ( |6.11 ) and so on. Continued 



iteration does not affect significantly the value we extract for the residue, but converges 
to an accurate measurement of the mass. The values in ( |6.13| ) were obtained in this way. 



With Zamolodchikov's S'-matrix the best fit of the data (with the same iterative 
procedure as before) gives the result 

A = 0.96 ±0.02; 

B = -1.10 ±0.02; (6.14) 
C = 1.14 ±0.02. 

The residue extracted from these data (1.29 ±0.01) is not consistent with that one of the 

\/3m 

amplitude a{f3). The situation does not improve even if we fix the coefficient of e 2~ to 
be that one predicted by Table 2, namely B = —0.158 and leave as free parameters for 
a best fit A and C. In this case, our best determination of A and C were A = 0.965 and 
C = —0.046. The curve is plotted in fig. 8 together with the data obtained from TCSA. 

7 Conclusions 

The S'-matrix proposed by Zamolodchikov for the subleading magnetic perturbation of 



TIM is a particular solution of the "Hard Square Lattice Gas" [^, fitted in such a way 
that it matches the physical picture coming from the TCSA data, i.e. the presence of 
two fundamental kinks and only one bound state thereof. But also the RSOS S'-matrix 
computed in sect. 3 is a particular set of Boltzmann weights for the "Hard Square Lattice 
Gas" and it reproduces the same features. The key difference is how the crossing sym- 
metry comes about: for the RSOS S-matrix, one finds a non-trivial charge conjugation 



matrix, eq. (|3.6| ), whereas for Zamolodchikov's proposal the crossing symmetry occurs in 



a standard way, eq. ( ^.5| ). Their asymptotic behaviour for large real values of (3 is also 



quite different: the RSOS 5- matrix goes to a definite limit while the S*- matrix proposed 
by Zamolodchikov is oscillating. For the RSOS case, the finite limit of the S'-matrix 



allows us to introduce generalized statistics of the kink excitations [O, na, 17 



The real discrimination between the two scattering theories proposed for the sublead- 
ing magnetic deformation of TIM is seen by comparing them with a numerical "experi- 
ment" , i.e. from the study of the finite-size corrections of the energy levels obtained by 
the Truncated Conformal Space Approach. We have investigated this problem in sect. 
6. The result suggests that the RSOS S'-matrix gives a more appropriate description for 
the scattering processes of the massive excitations of the model. Though, the interest- 
ing question as to what kind of system the S'-matrix of Zamolodchikov corresponds to, 
remains open. 
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